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This  technical  report  was  prepared  by  Harry  M.  Dobbins  and  Michael 
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number  7629/08/41.  It  is  the  final  report  for  the  period  January  1978 
to  October  1979  and  was  submitted  for  publication  in  December  1979. 
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SECTION  I 
INTRODUCTION 

1 .  STATEMENT  OF  THE  PROBLEM 

In  an  airborne  engagement  In  which  a  defending  platform  Is  subject 
to  multiple  simultaneous  attack,  decisions  must  be  made  concerning  which 
threats  must  be  countered  and  the  order  in  which  they  must  be  countered. 
These  decisions  can  be  based  on  the  time  of  arrival  of  the  Individual 
threats  and  the  probability  of  survival  of  the  defending  platform  based 
on  extrapolation  of  the  current  threat  trajectories. 

It  is  Important  to  note  that  the  errors  associated  with  extrapola¬ 
tion  of  threat  trajectories  are  a  function  of  two  major  error  sources. 

Measurement  uncertainties,  and 
Model  inaccuracies 

Given,  measurements  and  estimates  of  a  threats  position,  velocity, 
and  acceleration,  estimates  of  the  future  threat  state  can  be  made  by 
selecting  a  suitably  sophisticated  trajectory  model.  The  problem  lies 
In  selecting  a  model  which  is  useful  in  the  sense  that  it  Is  both 
computationally  efficient  and  reasonably  representative  of  the  process 
involved.  The  effective  use  of  this  extrapolation  approach  requires 
evaluation  of  the  error  Inherent  in  the  model.  In  using  this  approach 
the  estimate  of  the  error  must  be  Included  In  the  algorithm  Itself  or 
decisions  must  be  made  depending  on  the  specific  scenario  and  knowledge 
of  the  algorithm's  limitations. 

Estimation  of  error  places  additional  computational  burden  on  the 
system,  but  provides  a  more  generic  algorithm  In  the  sense  that  all 
scenarios  may  be  treated  by  the  same  algorithm.  The  second  approach 
requires  a  priori  knowledge  of  error  effects  for  classes  of  scenarios, 
the  algorithms  used,  and  the  usually  detailed  Identification  information. 
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Figure  1  Illustrates  a  typical  engagement.  The  predicted  miss 
distances  depend  on  the  assumptions  made  about  the  trajectory  and  noise 
In  measurements.  The  constant  velocity  assumption  (dashed  straight  line 
In  the  figure  from  missile  (1))  Indicates  the  trajectory  when  the  missile 
had  found  an  Intercept  course  and  Is  flying  at  essentially  constant 
velocity.  The  errors  Inherent  In  the  assumptions  would  not  provide 
accurate  miss  distance  estimates  for  missile  (1)  In  the  figure  which  Is 
not  on  an  Intercept  course  but  Is  still  correcting  to  an  Intercept 
course. 

The  constant  acceleration  assumption  (curved  line  In  the  figure 
from  missile  (1))  provides  an  Improved  estimate  of  miss  distance.  But 
now  noise  In  the  acceleration  estimates  and  actual  acceleration  due  to 
differences  from  a  constant  velocity  course  would  degrade  the  prediction 
for  missile  (2). 

Finally,  for  missile  (3)  both  assumptions  would  fall,  since  the 
threat  Is  not  on  an  Intercept  course  but  Is  executing  a  high  g  maneuver 
to  correct.  To  assume  a  constant  velocity  falls  to  take  Into  account 
the  high  normal  acceleration.  However,  to  assume  that  the  acceleration 
will  remain  constant  throughout  the  entire  trajectory  produces  an 
overcorrection.  The  assumptions,  error  sources,  and  possible  solutions 
to  these  problems  will  be  discussed  more  fully  in  Sections  2  and  3, 

Since  the  prediction  of  miss  distance  is  sensitive  to  both  the 
measurements  and  Inherent  assumptions  about  the  threat  trajectory,  a 
highly  sophisticated  threat  model  would  not  solve  the  problem.  Such  a 
model  would  not  only  Increase  the  computational  burden  of  the  system 
but  would  also  be  subject  to  Its  own  set  of  a  priori  assumptions  and 
measurement  noise. 
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An  alternative  solution  would  be  to  provide  an  estimate  of  the 
standard  deviation  for  a  number  of  successive  predictions.  If  the  threat 
trajectory  Is  accurately  determined  by  the  model,  the  differences  In 
prediction  will  be  relatively  small,  and  so  will  the  standard  deviation. 

If  the  acceleration  states  are  rapidly  changing  or  the  tlme-to-go 
value  Is  large  the  standard  deviation  will  be  correspondingly  large. 
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SECTION  II 
TECHNICAL  STUDIES 

1.  COMPARISON  OF  TRAJECTORY  ASSESSMENT  METHODS 

The  purpose  of  the  trajectory  assessment  function  is  to  evaluate 
and  rank  targets  with  respect  to  the  dangers  they  present  to  the  defend¬ 
ing  platform,  and  the  expected  time  of  intercept  or  closest  approach. 
These  functions  are  based  upon  the  state  estimates  of  position,  and  the 
velocity  and  acceleration  estimates  provided  by  the  sensors  and  tracking 
filters  which  can  be  used  to  predict  or  model  its  behavior. 

a.  Constant  Velocity  Assumption 

Time-to-go  and  predicted  miss  distance  are  based  on  the 
assumption  that  both  the  defending  platform  and  the  threat  missile 
maintain  a  constant  relative  velocity  vector.  The  predicted  miss  dis¬ 
tance  is,  then,  the  distance  at  the  point  of  closest  approach  and  the 
time-to-go  is  the  time  of  closest  approach.  This  calculation  is  some¬ 
what  simplified  by  adapting  an  intercept-plane  coordinate  frame.  This 
plane  Is  defined  by  the  relative  position  and  velocity  vectors  of  the 
defender-threat  pair.  The  aircraft  and  missile  velocities  normal  to  the 
intercept  are  equal,  by  definition.  Furthermore,  If  all  significant 
accelerations  occur  In  the  intercept  plane  the  orientation  of  the 
intercept  plane  will  not  change  as  the  scenario  unfolds. 


Figure  2  illustrates  the  geometric  relationships  involved, 
when  It  is  assumed  that  the  defending  platform  and  threat  maintain  a 
constant  velocity.  The  velocity  of  the  threat  relative  to  the  platform 
Is 


V 


t 


+  V  ) 


2  ,!/2 


(1) 
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where  r  is  the  velocity  along  the  line  of  sight  (LOS)  and  Vp  Is  the 
velocity  normal  to  the  line  of  sight  and  is  given  by 

Vn  =  r(<)>2  +  e  cos  (j*)1^2  (2) 

The  heading  error  (a)  Is  defined  to  be  the  angle  of  the  LOS  and  the 
relative  velocity  vector  as 


sin  a  =  Vn/Vt  (3) 

From  the  geometry  of  the  figure, 

sin  a  =  XQ/r  (4) 

where  XQ  is  the  predicted  miss  distance,  or  point  of  closest  approach, 
and  r  is  the  range.  Then  n  can  be  calculated  from 


X 


o 


=  r  V/Vt 
n  t 


and  the  time-to-go  is,  then. 


*  2 

-r  r/v/ 


(5) 


(6) 


In  practice,  however,  a  guiding  threat  does  not  maintain 
a  constant  velocity  relative  to  the  LOS.  Small  deviations  from  the 
true  collision  course  would  induce  changes  in  V  and  r,  due  to  the 
missile  guidance,  even  if  the  total  velocity  did  remain  constant.  Since 
the  total  velocity  actually  changes  due  to  boost  and  drag,  the  intercept 
course  does  change. 
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Assuming  these  changes  in  velocity  are  relatively  small  and 
the  greatest  source  of  error  will  be  the  measurement  or  estimate  error, 
then,  from  Equation  5, 


AXft  AV  A„ 
_ o_  _ _ n  ,  _r 

Xft  "  Vn  r 

o  n 


and,  from  Equation  6, 


AT  A 
go  _  r 


Then,  for  a  small  heading  error  a,  V  will  be  large,  compared  to  aV. 
and  r,  and  r  will  be  large  compared  to  ar  and  ar.  However,  the  normal 
velocity  Vn  will  be  of  the  same  order  of  magnitude  as  aVn>  By  definition 
of  constant  velocity  intercept  trajectory,  Vp  =  0.  Small  errors  in  the 
estimates  of  Vn  will  dramatically  affect  the  miss  distance  prediction. 


b.  Constant  Acceleration  Assumption 

The  approach  taken  and  discussed  by  the  Charles  Stark  Draper 
Laboratories  (Reference  2)  includes  the  use  of  constant  velocity  and 
assumptions  of  both  small  axial  and  normal  acceleration.  Then  the  miss 
distance  and  time-to-go  are  expressed  as 


■  Vgo 


Vgo 


Tnrt  *  (r  •  VJ 
go  _o — a 

iCl 
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where 

rQ  Is  the  missile  initial  position 

VQ  is  the  Initial  velocity 

A  is  the  initial  and  assumed  constant 
acceleration 

Equations  9  and  10  will  provide  exact  results  in  the  event  that  a 
coasting  missile  approaches  the  defending  platform  along  its  implied 
constant  velocity  vector.  Any  deviation  from  this  set  of  circumstances 
would  result  in  an  undetermined  error  associated  with  the  estimator. 
Reference  2  details  a  more  complicated  version  of  Equation  9  which  would 
provide  an  enhanced  computation  of  multi-aspect  threat  arrivals,  but 
drag  and  staging  will  still  result  in  accelerations  not  modeled  by  the 
estimator  (Reference  3).  This  approach  also  suffers  because  it  does 
not  attempt  to  model  the  missile  guidance.  The  normal  acceleration  will 
not  remain  constant,  due  to  guidance,  but  will  tend  toward  zero  as  the 
intercept  course  is  attained. 

c.  Decreasing  Acceleration  Assumption 

To  account  for  decreasing  normal  acceleration  as  the  intercept 

course  is  attained  the  following  assumptions  are  made.  The  normal 

acceleration  decreases  linearly  with  time  and  at  T  the  normal  velocity 

9® 

reaches  zero.  Then  the  normal  velocity  becomes 

Vn  -/V'-V  t  <"> 

Vn  '  V  ‘  Hi  t2't5o  *  Vo  <12> 

and  the  miss  distance  becomes 

vn 

h  -f  \  \  <”> 

v_ 
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Then  the  tlme-to-go  can  be  calculated  from  the  radial  acceleration  and 
velocities  (assumed  constant)  by 


Although  this  method  attempts  to  account  for  missile  guidance,  no 
specific  guidance  Is  modeled,  and  the  effects  of  staging  and  nonlinear 
drag  effects  are  still  not  modeled. 

2.  COMPARISON  OF  METHODS 

All  of  the  methods  discussed  are  based  on  a  polynomial  approxima¬ 
tion  of  the  threat  trajectory  and  extrapolation  to  the  point  of  closest 
approach.  The  constant-velocity  method  is  based  on  linear  or  straight- 
line  approximation,  and  does  not  account  for  acceleration.  The  constant- 
acceleration  method  uses  a  quadratic  approximation,  but  does  not  account 
for  guidance.  The  decreasing  acceleration  method  Is  based  on  a  cubic 
approximation,  which  Is  reduced  to  a  quadratic  expression  by  making 
assumptions  about  the  threat  trajectory. 
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The  major  problem  associated  with  each  of  these  methods  Is  a 
divergence  of  the  predicted  miss  distance  from  the  actual  miss  distance 
for  large  tlme-to-go.  Since  the  value  of  function  of  t  and  Its  derivl- 
ties  are  specified,  the  error  associated  with  this  type  of  extrapolation 
has  an  error  of  the  form  (Reference  6). 


e  =  F(t)  -  P(t) 


(18) 


e  «  F(t)  - 


AoT"-1 


.+  An)  (19) 


If  P(t)  Is  an  accurate  representation  of  F(t),  the  error  grows 
In  accordance  with  the  power  of  the  number  of  steps  extrapolated  forward 
multiplied  by  the  error  In  the  coefficients.  Figure  2  illustrates  the 
error  associated  with  the  methods  discussed.  As  indicated,  the  constant 
velocity  assumption  generally  gives  poorer  estimates  than  either  a  con¬ 
stant  or  a  decreasing  acceleration  assumption.  The  constant  acceleration 
assumption  provides  the  best  estimate  when  the  threat  must  actually  main¬ 
tain  constant  acceleration  to  arrive  at  the  point  of  closest  approach, 
otherwise  the  decreasing  acceleration  assumption  provides  a  better  esti¬ 
mate.  The  major  disadvantage  to  using  these  approaches  Is  the  lack  of 

inherent  measure  of  validity.  That  Is,  there  Is  no  measure  of  at  and 

go 

a  (Reference  3).  However,  since  the  miss  distance  must  be  repetitively 
evaluated  In  real  time,  the  average  and  standard  deviations  associated 
with  the  miss  distance  can  also  be  calculated  by 


where  X0^  Is  the  1th  estimate  of  miss  distance 

X0  Is  the  average  value  of  the  last  N  estimates 
av 

N  Is  the  number  of  estimates 

o  Is  the  standard  deviation  over  the  last  N  estimates 
xo 


(20) 
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Actual 

Constant  A„ 
n 


Predicted 
Constant  An 

Decreasing  An 
Constant  Vn 


Figure  3.  A  Comparison  of  Actual  Constant  An  Missile  Trajectory  and 
Predicted  Miss  Distances  With  Constant  Velocity  Normal  to 
LOS,  a  Constant  A  ,  and  Decreasing  An  In  Alrcraft-to-MIssIle 
LOS  Frame 


AFAL-TR-79-1 242 


The  standard  deviation  calculated,  not  only  provides  a  measure  of 
validity  of  the  model  used,  but  can  also  be  used  directly  in  the  calcu¬ 
lation  of  probability  of  survival. 

3.  PROBABILITY  OF  SURVIVAL  WHEN  GIVEN  NO  ACTION  (PS/NA) 

When  given  that  a  missile  detonates  at  a  distance  x  from  an  air¬ 
craft,  we  can  assum,e  the  probability  of  survival  depends  on  the  detona¬ 
tion  distance  x  and  the  warhead  lethality  range  rQ  and  express  it  as 

PS/NA  =  1  -  e"x2/ro2  (21) 

For  simplicity,  assume  that  r0  is  not  a  function  of  the  angular 
position  of  the  missile  relative  to  the  aircraft's  velocity  vector. 

Then  the  value  of  rQ  defines  the  1  sigma  radius  for  a  given  lethality. 
Without  identification  information  the  value  of  r0  would  be  unknown. 

In  this  case  r0  would  be  selected  as  a  maximum  desirable  keep-out-of¬ 
range, 

A  defensive  system  does  not  have  absolute  knowledge  of  what  the 
value  of  x  is  before  the  detonation  of  the  missile.  The  defensive  sys¬ 
tem  can  make  an  estimate  of  the  point  of  closest  approach  and  can  make 
an  estimate  of  the  time-to-go  before  the  missile  is  at  the  point  of 
closest  approach.  The  error  between  the  estimated  point  of  closest 
approach  and  the  actual  point  of  detonation  results  because  the  defensive 
system  cannot  accurately  predict  the  dynamic  motion  of  the  missile  for 
long  periods  of  time.  In  general,  the  higher  the  value  of  the  standard 
deviation  of  x  before  the  missile  Is  at  the  point  of  closest  approach, 
the  more  uncertain  the  estimate  of  the  point  of  closest  approach  (xQ) 
will  be. 

Assume  the  probability  density  of  a  missile  having  a  detonation 
•iiscance  of  x  to  be:  /x“ x0\  2 

dP(x)  .  'V^T/  (22) 

dx  -we 
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PREDICTED  MISS  DISTANCES 


Figure  4.  Probability  of  Survival  as  a  Function  of  Miss  Distance  for 
Various  Lethality  Ranges 
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where  ox  is  the  standard  deviation  of  the  value  of  xQ,  and  N  is  the 

normalization  coefficient.  If  a  =  0,  then 

*o 


The  value  of  N  can  be  determined  because  the  total  probability 
over  all  values  of  x  from  0  to  °°  must  be  one 


(23) 


ot 

/ 


Let 


Then 


and 


1  =  N  /  e 
o 

x-x_ 


-(£*) 

o'. 

dx 


dx  =  __o  dy 
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xo 
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(24) 
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The  error  function  Is  defined  as 


ERF  (*?)-  ? 

f  e‘^2  dy 

(29) 

0 

ERF  («)  *  1 

(30) 

1  =  Ct  TT 

N  —22 - 

2  \ 

iltERF  (£.)} 

*0 

(31) 

N  =  2 

\  ,  {  i 

+  ERF(  ox°  )} 

(32) 

The  probability  of  survival  (given  that  no  action  Is  taken)  Is, 
then,  the  Integrated  product  of  the  probability  density  of  the  point  of 

detonation  occurring  at  x  and  the  probability  of  surviving  the  detonation 
at  x,  and  can  be  written 


PROBABILITY  OF  SURVIVAL 
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Let  x„  and  a  be  expressed  In  units  of  r  ,  so  that 
ox  o 


Figures  5  and  6  illustrate  the  effect  of  uncertainty  or  model  error 
for  various  miss  distances.  Figure  5  shows  that  as  estimated  miss  dis¬ 
tance  increases  the  probability  of  survival  also  increases.  Figure  6 
shows  that  as  the  estimate  Improves  (the  modeled  behavior  of  the  threat 
resembles  the  actual  behavior),  better  estimates  of  probability  of 

survival  are  generated.  As  ox  becomes  very  large  compared  to  miss 

0 

distance  the  estimated  probability  of  survival  becomes  PS/NA  =  1-e  . 

Figure  7  gives  a  contour  plot  of  the  probability  of  survival  for 

miss  distance  and  cx  .  Note  that  as  the  predicted  miss  distance  and 

o 

ax  increase,  the  probability  of  survival  estimate  is  less  sensitive  to 
o 

variations  In  miss  distance. 

4.  THREAT  PRIORITIZATION 

The  priority  of  trie  threat  is  based  on  the  probability  of  survival 
and  Tact  and  the  time  or  range  at  which  a  defensive  action  can  be  taken. 

The  time  at  which  defensive  action  can  begin  can  be  calculated 

from 

Tact  '  <K  -  *act)A 
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For  the  constant  velocity  assumption,  where 
R  Is  the  range, 

Ract  the  range  at  which  action  can  begin,  and 
VL  Is  the  velocity  along  the  line  of  sight. 

Or,  In  the  constant  acceleration  case 


where 

Is  the  acceleration  along  the  LOS. 

Then  the  priority  can  be  expressed  as 

P  -  (1  -  PS/Nft)  ERF  ((Tgo  -  Tact)/Tg0) 

The  error  function  rapidly  Increases  towards  unity  as  TgQ  becomes 
less  than  Tact,  In  which  case  the  priority  is  then  determined  by  the 
probability  of  kill  due  to  the  threat. 
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SECTION  III 
CONCLUSION 

The  methods  developed  for  deterministic  extrapolation  of  threat 
trajectories  to  perform  probability  of  survival  and  prioritization  have 
several  positive  features: 

-  relative  simplicity 

-  Indication  of  relatively  good  accuracy  over  a  range  of  test 
scenario  (References  1  and  2). 

The  problems  associated  with  use  of  these  methods  are  error  sources: 

-  measurement  errors  and  estimation  inaccuracies 

-  minor  model  inaccuracies  (the  missile  guidance  laws  are  not 
directly  modeled)  which  cause  divergence  of  the  extrapolated  trajectory 
from  the  actual  trajectory. 

-  major  model  inaccuracies,  such  as  staging  and  nonlinear  drag 
effects. 

The  major  pro^em  associated  with  deterministic  extrapolation  has 
been  the  lack  of  statistical  measurement  of  the  validity  of  the  model. 

The  statistical  information,  standard  deviation  of  miss  distance,  can  be 
directly  calculated  from  N  sequential  estimations  and  directly  used  to 
determine  probability  of  survival  by  the  method  developed  in  Section  II. 3. 
This  reduces  the  sensitivity  of  the  probability  of  survival  calculation 
when  the  trajectory  model  is  inaccurate  or  is  giving  rapidly  N  changing 
results. 

The  model  inaccuracies  associated  with  missile  guidance,  staging 
and  drag,  would  require  detailed  and  complex  models.  This  type  of  model 
would  either  require  some  type  of  identification  information  or  methods 
of  determining  the  missile's  guidance  gain  in  real  time.  There  are 
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nonlinearities  associated  with  staging  and  uncertainties  in  predicting 
when  it  would  occur.  Since  the  current  techniques  provide  reasonable 
results,  the  additional  complexities  and  Increased  computational  burden 
of  a  more  detailed  model  would  seem  to  preclude  the  possibility  of  a 
real-time  solution  using  a  more  complex  model  at  this  time. 

The  measurement  errors  and  state  estimation  inaccuracies  would  be 
Inherent  in  any  system,  but  they  become  more  critical  when  polynomial 
extrapolation  techniques  are  used,  since  the  error  in  estimation  grows 
as  the  product  of  the  error  and  a  power  of  the  number  of  time  steps 
extrapolated  forward. 

Methods  for  reducing  the  effect  of  this  error  could  include  various 
data  smoothing  techniques,  such  as,  the  least-squares  polynomial  fit  or 
the  Chebyshev  min-max  polynomial  fit.  These  techniques  are  attractive 
since  the  number  of  data  points  necessary  to  calculate  the  standard 
deviation  required  for  the  probability  of  survival  estimation  could  also 
be  used  for  data  smoothing.  This  would  reduce  the  more  critical  errors 
associated  with  the  polynomial  approximation. 

Few  data  points  are  required  since  the  degree  of  the  polynomial 
approximation  is  low.  For  example,  only  four  data  points  would  be 
required  to  generate  a  Chebyshev  min-max  parabola,  which  is  the  highest 
degree  polynomial  discussed. 

Although  these  techniques  reduce  the  measurement  error,  they  do 
add  computational  burden  and  do  not  eliminate  model  inaccuracies.  Since 
model  inaccuracies  exist,  repetitive  calculations  are  necessary.  Addi¬ 
tional  trade  studies  would  be  required  to  determine  the  benefit  associated 
with  Increased  accuracy  versus  Increased  computational  burden  for  a  parti¬ 
cular  system. 
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APPENDIX 

Figures  Illustrating  Nominal  Test  Cases 
Using  the  Algorithms  Developed  Using  the 
Constant  Velocity  Assumption  and  Acceleration 
Assumptions 
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Figure  7.  Nominal  Trajectory  No.  1,  Head-On  AC  at  45°  to 
Missile  Heading 
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Figure  9.  Tlme-to-Go  for  Trajectory  No.  1 
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Figure  10.  Time-to-Go  for  Trajectory  No.  2 
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Figure  13.  Probability  of  Survival  for  Trajectory  No.  1 
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Figure  15.  Priority  Ranking  for  Trajectory  No.  1 
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